Vortex lattices for ultracold bosonic atoms in a non-Abelian gauge potential 
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The use of coherent optical dressing of atomic levels allows the coupling of ultracold atoms to 
effective gauge fields. These can be used to generate eff'ective magnetic fields, and have the potential 
to generate non-Abelian gauge fields. We consider a model of a gas of bosonic atoms coupled to a 
gauge field with f/(2) symmetry, and with constant effective magnetic field. We include the effects 
of weak contact interactions by applying Gross-Pitaevskii mean-field theory. We study the effects of 
a f/(2) non-Abelian gauge field on the vortex lattice phase induced by a uniform effective magnetic 
field, generated by an Abelian gauge field or, equivalently, by rotation of the gas. We show that, 
with increasing non-Abelian gauge field, the nature of the groundstate changes dramatically, with 
structural changes of the vortex lattice. We show that the effect of the non-Abelian gauge field is 
equivalent to the introduction of effective interactions with non-zero range. We also comment on the 
consequences of the non-Abelian gauge field for strongly correlated fractional quantum Hall states. 

PACS numbers: 03.75.Lm, 73.43.Nq, 03.75.Kk 



I. INTRODUCTION 

Atomic Bose-Einstein Condensates (BECs) offer the 
possibility to study the physics of quantised vortex lines 
with unprecedented precision and control [TJ [5] . Exper- 
iments on rapidly rotating gases [3HS] have allowed de- 
tailed studies of remarkable states such as large arrays 
of vortices, forming vortex lattices. Their static features, 
dynamics and response to periodic lattice potentials have 
been investigated. 

As an alternative to rotation, one can use the dress- 
ing by coherent optical fields to create an effective C/(l) 
gauge potential which simulates the orbital effects of a 
magnetic field on a charged particle [?]• Using such op- 
tically induced gauge fields, the formation of quantized 
vortices in a rubidium condensate has been demonstrated 
in pioneering experimental work |5]. Optically induced 
gauge potentials are not limited to Abelian gauge fields, 
but can be naturally extended to the non-Abelian case 
[7]. There exists a variety of proposed ways to gener- 
ate non-Abelian gauge potentials, both in the continuum 
[9, lOj and in lattice-based systems [ITj 

In this paper, we study the consequences of a U(2) 
non-Abelian gauge field on the groundstate of a weakly 
interacting atomic BEC. We focus on a gauge-field con- 
figuration in which the effective magnetic field is constant 
in space, and for which there exists a simple exact solu- 
tion for the single particle wavefunctions [HHH]- In the 
case of a uniform Abelian magnetic field (or for uniform 
rotation of the gas) the spectrum has the Landau level 
structure. In this case, for weak repulsive interparticle 
interactions the bosonic atoms occupy the lowest energy 
Landau level and the mean-field ground state is well- 
known to be a lattice of vortices with triangular symme- 
try. It is interesting to ask how the addition of a constant 
non-Abelian magnetic field affects this groundstate. 

We show that the effects of the additional non-Abelian 
magnetic field can be understood in terms of a change 



in the effective interatomic interaction potential, and 
are equivalent to the effects of an interaction potential 
with non-zero range. Specifically, the effects of the non- 
Abelian gauge field are completely encoded on the Hal- 
dane pseudopotentials that describe the interatomic in- 
teractions in the lowest energy single-particle states. We 
study the consequences on the condensed vortex lattice 
phases by applying Gross-Pitaevskii mean-field theory to 
the weakly interacting gas. We show that the nature of 
the ground state changes dramatically, with structural 
changes in the symmetry of the vortex lattice brought 
about with increasing non-Abelian gauge field. We show 
that these changes are precisely analogous to the intro- 
duction of a long-range interaction, as has previously 
been studied in the context of dipolar interactions [16j . 

The paper is organised as follows. In Sec. |TT]we intro- 
duce the non-Abelian gauge field. In Sec. |III| we introduce 
the interaction Hamiltonian and the corresponding Hal- 
dane pseudopotentials. In Sec. |IV| we present the vortex 
lattices. In Sec. [V] we comment on fractional quantum 
Hall states in the system. Sec. |VI| contains our conclud- 
ing remarks. Some calculational details are relegated to 
Appendices. 



II. A NON-ABELIAN GAUGE FIELD 

The Hamiltonian for a non-relativistic particle of mass 
m and charge q reads 
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where p is the particle momentum and A is a vector po- 
tential which produces the magnetic field B ~ V x A. 
We are interested in the case of a non-Abelian vector 
potential which is written in the form 



A — AxUx + Ayiiy -|- AzUz 



(2) 



and the components A^ are matrices acting on a set of 
states which, in the present context, correspond to a set 
of degenerate dressed states [7]. In the simplest case there 
are two such internal states, which we label by u =t, J,, 
and Ay, arc 2x2 Hermitian matrices. The gauge group is 
then t/(2) = C/(l) x SU{2) and thus contains the standard 
f/(l) potential and a SU{2) part. The corresponding 
magnetic field is [T7] 



B = W X A-i^Ax A. 
n 



(3) 




This relation shows that, in the case that the compo- 
nents of A are non-commuting, a nonzero magnetic field 
is obtained even from a uniform vector potential. 

Let us consider a uniform magnetic field perpendicular 
to the plane B = BzUz, where i?^ is a 2 x 2 hermitian 
matrix. The magnetic field can be assumed to be di- 
agonal by an appropriate choice of basis and we write 
Bz = B{X + 20^(Jz), where X is the identity matrix, Cz is 
the diagonal Pauli matrix, and /3 is a parameter control- 
ling the size of the non-Abelian part of the field. The first 
Abelian term in such a magnetic field can be produced 
in the standard way by the rotation term in Eq. (Is]). 
The second term of the magnetic field is produced by a 
vector potential whose components are non-commuting 
constant matrices, through the second term on the right- 
hand side of Eqn. (Sl). The complete non-Abelian vector 
potential may be chosen in the form 



A 




(4) 



where ax,Uy are Pauli matrices. The constants in the 
magnetic field and the vector potential are related by /3 = 
f3' /{£bB), where we have introduced the magnetic length 
£b = \/fi/{qB). We have chosen the Landau gauge for 
the Abelian part, while the second, non-Abelian term 
has been chosen to stand in analogy to the symmetric 
gauge. Up to an overall gauge transformation this gauge 
potential is equivalent to the gauge potentials studied 
in Ref. [T2J [T3] and for the symmetric case \a\ — \b\ in 
Ref. fT5] . A method for implementing such a non-Abelian 
gauge field is described in Ref. [13[. 

We introduce the standard creation and annihilation 
operators a^ , a for the Landau level problem and the spin 
ladder operators a± — (cFx i ay)/2, and then the Hamil- 
tonian is written in the form 



H = huJr 



i'^a + V2(3{a^a+ +ao-_) + - +/3^ 



, (5) 



where Wc = qB/m is the cyclotron frequency. This is 
equivalent to the Jaynes-Cummings Hamiltonian,I8j and 



FIG. 1: The energy levels m for n = 0, 1, 2, 3 as functions of 
j3. We find that Ei is the minimum energy for < /3 < -\/3 
while E2 is the minimum for \/3 < P < \/5- 



its spectrum is 



Eq = hwc 
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(n±^+pA, /3„ = V8/32n+I, n > I. 



(6) 



We are interested in the ground state as a function of 
the parameter /3. Since E^ > E^ we will only discuss 
En = E~ . Fig. [I] shows the first few energy levels as a 
function of /?. One finds that Ei is the lowest energy for 
< /3 < -\/3, while the ground state energy is En for 
i/2n- I < /3 < y/2n + 1 when n > 1. 

Let us denote the Landau level states by (fin.k where n 
is the Landau level index and k is an index for the de- 
generate states in each Landau level. The two spin states 
will be denoted by the symbols t and l so the product 
states of the Landau level states with the spin states are 
denoted as 4>n,k,^, 4>n,k,i- The normalized eigenfunctions 
for the energies En are 
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In Appendix A, we give explicit expressions for the spa- 
tial dependence of the wavefunctions 4>n,k in the periodic 
geometry used below for our numerical calculations. 



III. INTERACTION HAMILTONIAN 

We assume that the bosons are interacting. We use a 
second-quantized description, and write t/iJ. (r) and ■00- {r) 
for the creation and annihilation operators of a particle 
with spin a at position r. The interaction Hamiltonian 
for two-body interactions has the form 



Ht = 



\ E [9..'{\r-^\) ^l{7'Hl{r)Mr^i^^'ir') d'rd'r^ 



(8) 



We have assumed that the interaction potential depends 
on the distance between the particles and it may also 
depend on the spins. We have made the additional as- 
sumption that the spins of the bosons do not change due 
to scattering. We suppose a contact interaction potential 
of the form g^^,[\f— r'\) = gS{f— f) and obtain 



H, 
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V'^/V'ffV'tTV'CT' d?r. 



(9) 



(It is straightforward to retain the spin-dependence to 
the interactions, which will introduce additional param- 
eters to the model. We restrict attention to the spin- 
independent case to simplify presentation.) 

We consider, first, the effect of interactions within the 
basis of single-particle eigenstates of definite angular mo- 
mentum y\. We assume that the interactions are small 
compared to the energy level spacings between states 
with different quantum numbers n, so we focus only on 
the scattering of pairs of particles with the same value of 
n, but between different angular momentum states: that 
is (n, mi), (71,7712) — )■ (77,7711), (77,7772). Owing to the ro- 
tational symmetry of the (contact) interactions, all two- 
body scattering processes ^ preserve the total angular 
momentum of the pair of particles, m'l -\- m'2 = mi + 7772. 
As a result, the two-particle states with definite rela- 
tive angular momentum, 777, are eigenstates of the in- 
teraction Hamiltonian; their eigenvalues are the Haldane 
pseudoDOtentials|191 I20j. Vm \ which encode the entire 
properties of the interactions within the energy band la- 
belled by 77. For the states ^ these take the form 
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(10) 



We use here and in the following as a unit of length the 
magnetic length Ib- 

For bosons, symmetry of the two-particle wave func- 
tions means that only even 771 = 0, 2, 4, . . . can contribute. 
For the usual case of rapidly rotating bosons in the low- 
est Landau level, n = 0, interacting via contact interac- 
tions [2], the only non-zero pseudopotential is Vq. This 
represents the situation in which interactions have the 
shortest possible range. Explicitly introducing a long- 
range interaction potential, such as the dipolar potential 
with g{r) ex; l/r^ at large distances, leads to non-zero val- 
ues of pseudopotentials with m > 0[TB]. We find, from 
an analysis of (10 1, that when the non-Abelian gauge 
field is present, /? ^ 0, a set of pseudopotentials with 
771 > are also non-zero even for contact interactions 
g{\r\) oc 5{r). We therefore argue that the effect of the 
non-Abelian gauge field on the properties of the lowest 
energy state is equivalent to the introduction of effective 
interactions with non-zero range. 

In the case of dipolar interactions, these changes to 
the effective interaction were found to lead to changes in 
the lowest energy vortex lattice phase [TB]. In the fol- 
lowing we investigate the effects on the vortex lattices 



due to the changes in effective interaction caused by the 
non-Abelian gauge field. We first give results of numeri- 
cal simulations and we then make the connection to the 
calculated values of the pseudopotentials Vm ■ 



IV. VORTEX LATTICE PHASES 

To study vortex lattices of an infinite system, it is con- 
venient to work in a rectangular geometry with periodic 
boundary conditions in the two spatial co-ordinates, that 
is with the topology of a torus. We consider a rectangular 
cell of size L^ and Ly in the x and y directions respec- 
tively. This spatial periodicity imposes the condition 



i^iy = 27r47Vv, 



(11) 



where iVy is a positive integer equal to the number of 
states in the cell for each Landau level. This can be in- 
terpreted as the number of vortices within the cell. De- 
pending on the symmetry of the vortex lattice phase (de- 
scribed within mean field theory below), the aspect ratio 
of the cell Ly/L^ must be chosen to match the natural 
periods of the vortex lattice. 

For small enough values of the interaction strength 
gp <C hujc, where p = N/L^Ly is the mean particle den- 
sity, we can assume that the interaction potential does 
not mix excited states of the non-interacting Hamilto- 
nian. Therefore the ground state of the interacting sys- 
tem is found by minimizing the interaction energy within 
the space of the wave functions (l7| for a specific n. 

We make the ansatz 



*n = y^^Ckip„,k 



(12) 



where c^ are complex coefficients, and ■(/'„, fc are the single 
particle wavefunctions (l7| for the periodic geometry (see 
Appendix A). The choice for the Landau level index 77 
depends on the value of /3. We minimize the energy 



Er = 



E y^ 



fci,fc2,/C3,fc4 C/jj^ Cj,^ Cfe3 Cfe^ 



(13) 



^1,^2,^3.^4 



(the form of the interaction potentials Vk-^^k2.ka,k4, is given 
in Appendix IbI) in the subspace of wave functions ( 12 ) 
for a given number of particles 



^ = E 



CkCk. 



(14) 



The mean density is then p = -^ ^p , and the condition 
that mean-field theory provides an accurate description 
is that N/Nv > 1 0. 



We find minima of the energy ( 13 1 numerically us- 



ing the Fletcher-Reeves-Polak-Ribiere method ||22j which 
uses derivatives of the function to be minimized. The 



number of basis wave functions in the ansatz ( 12 ) is Ny 
and we used values up to iVy = 64 for which the method 





FIG. 2: Contour plots for the particle density of a triangular 
vortex lattice for j3 = 1.0 (aspect ratio a^ — \/3). We present 
one cell which contains AV = 2 vortices, (a) The total particle 
density (in units of the mean particle density) . The particle 
density is nowhere zero, since the single particle wavefunctions 
involve contributions from more than one Landau level, as 
shown in Eqs. (7|, ( |12[ ). (b) The density for particles in the 
lowest Landau level (in units of the mean particle density in 
the lowest Landau level), (c) The density for particles in the 
first Landau level (in units of the mean particle density in the 
first Landau level). 



rapidly converges to a minimum. Setting iVy to a certain 
value we fix the number of Landau level states which are 
introduced in the area [L^:, 0] x [0, Ly], and the cell area is 



fixed through Eq. ( 11 ). As a second parameter we choose 



the aspect ratio a^ = Ly/L^ of the rectangular cell. We 
have explored the minima of the interaction energy in 
the space of the two parameters iVy, a^. For every state 
under consideration we evaluate the quantity 



^=|^(4'^^v), 



(15) 



where Ej is the interaction energy and N is the number 
of particles in a cell with Ny vortices, from which the 
chemical potential is /i = ^ = &p. 

For /3 = the ground state wave functions are the 
lowest Landau level states (/>o,fc,^, and we can verify that a 
triangular vortex lattice (identical to Fig. la in ReL [23 ) 
is the lowest energy state with h — 1.596, as expected by 
the equivalent problem in the Ginsburg-Landau model 
for type-II superconductors [MJ [2S] . 

For < /3 < -\/3 the ground state wave functions are 
obtained by setting n = 1 in Eq. (It]), and they are su- 
perpositions of the lowest and first Landau level states. 
Minimizing ( 13 ) we find a triangular lattice as the lowest 



energy wave function for the range < /3 < 1.13. Note, 
however, that the groundstate wave function depends on 
/3 explicitly through the coefficients of 0o,fc,i and <j)i^k,-\ in 
Eq. (It]). In Fig. [2] we present plots of the particle density 
for P = 1.0. We plot the total density as well as the den- 
sities for the two Landau levels of the wavefunction. As /? 
increases the contribution of 4>i^k,t in the wave function 
increases as indicated in Eq. ^ . The contribution to the 
particle density from the first Landau level wave function 
^i,fe,t has maximum values where the vortex centers for 
4'o,k,-\ are located. These remarks suffice in order to fol- 
low the change of the particle density in the lattice as 




FIG. 3: Contour plots for the particle density of a square 
vortex lattice for /3 = 1.5 (aspect ratio Or = 1). We present 
one cell which contains AV = 1 vortex, (a) The total particle 
density. The particle density is nowhere zero, since the single 
particle wavefunctions involve contributions from more than 
one Landau level, (b) The density for particles in the lowest 
Landau level, (c) The density for particles in the first Landau 
level. 



/? increases. Note that we have density minima, but no 
zeros of the density, in Fig. [2^. 

For 1.13 < /3 < -\/3 our numerical method converges to 
a square lattice which has lower energy than any other 
state we investigated. Fig. |3] shows the result of the en- 
ergy minimization in a single unit cell for a square lattice 
for P = 1.5. 

A transition from a triangular to square and other vor- 
tex lattices has been investigated for the case of long- 
range dipolar interactions in the LLL |16j . In order to 
quantify the contribution of non-local interactions we cal- 
culate the values of pseudopotentials Vm ■ It is conve- 
nient to define as a control parameter the ratio of the two 
first nonzero pseudopotentials at a certain Landau level 



V2 



(16) 



.(1) 



We have simplified the notation setting Vm = Vm' ■ If 
Vm — for TO > 2 a transition from triangular to square 
lattice occurs at a = ai = 0.0865 [53]. For rt = 1 in 
the present model we find that Vm — for to > 2, and 
the value a = ai is obtained for j3 — 1.13. Therefore, 
the transition to a square lattice for (3 > 1.13 reported 
here is in perfect agreement with the results of Ref. [33] . 
The parameter a increases with /3, and for /3 = v3 we 
have a — 1.118. In further agreement of the present 
results to the results of Ref. [23] , no transition to a vortex 
phase different than the square lattice is expected for 
these parameter values. 

When (3 increases to values -s/S < (3 < \/5 the ground 
state of the non-interacting problem is a superposition 
of the first and second Landau levels obtained from wave 
function (l7| with n = 2. Fig. [4] shows the particle density 
of the ground state for (3 ~ 2.0. The picture is similar for 
the whole range of values of f3. The particle density has 
vortex minima arranged in a triangular lattice and these 
are surrounded by density dips. The clusters of vortex 
and density dips may be called bubbles and we thus have 
a bubble crystal phase in Fig. [2] 

Following Refs. [TOl [23 we may characterize the vortex 
lattice in Fig. [4] as a g = 4 bubble as each unit cell con- 





FIG. 4: Contour plots for the particle density oi a. q — 4 
bubble crystal lattice for /3 = 2.0 (aspect ratio a^ — \/3). 
We present a cell which contains Ny = 8 vortices, (a) The 
total particle density. The particle density is zero on the sites 
of a triangular lattice, (b) The density for particles in the 
first Landau level, (c) The density for particles in the second 
Landau level. 



tains 4 vortices. The particular bubble phase of Fig. |4] 

is similar to that in Fig. le in Ref. ^3]. The Haldane 

pseudopotentials for /3 = \/3 are a = V2/Vb — 0.65 

and V4/V0 = 0.16, where we use the simplified notation 

(2) 
Vm — Vm ■ The vortex lattice for this value of a gives 

a q = 4 bubble state in the case of dipolar interactions 
[IB] as well as in the model of Ref. [23] • The present 
results are thus in agreement with the preceding reports 
showing that the non-Abelian gauge potential has an ef- 
fect analogous to long-range interactions on the vortex 
lattice phases. 

Note that we find no evidence of the appearance of 
"stripe crystal" phases, which appear for 0.12 < a < 
0.70 for the model of Ref. [23]. For the present model 
we have shown that < a < 0.12 for < /3 < -\/3 
and a > 0.70 for /3 > \/3 as a result of the transition 
for the groundstate in the noninteracting model at /3 = 
v3. A direct transition from square lattice to bubble 
crystal at /? = VS is therefore consistent with the results 
of Ref. jm. 



V. CONSEQUENCES FOR FRACTIONAL 
QUANTUM HALL STATES 

The above mean-field studies are valid in the regime of 
high filling factor i^ = N/Ny ^ 1 [2]. For small values 
of the filling factor, the vortex lattice phases can be re- 
placed by strongly correlated phases, which are bosonic 
analogues of fractional quantum Hall (FQH) states. For 
atoms in a uniform Abelian gauge field at z/ = 1/2 
and interacting with contact interactions, the ground- 
state at ly — 1/2 is the Laughlin state[5S|, and those at 
1/ = 1,3/2,2 are well described by the Moore-Read and 
Read-Rezayi statesI27j. Variations in the (ratios of the) 
Haldane pseudopotentials can lead to changes in the na- 
ture of the groundstate. 

Our result, Eqn (10), provides the values of the pseu- 



dopotentials for the non-Abelian gauge potential studied 
here. By combining these results with existing numerical 



studies of the effects of variations of the pseudopotentials 
on strongly correlated phases [Ml [23l [28l [29] we can de- 
duce consequences of the non-Abelian gauge field for the 
FQH states. 

For filling factor v = 1/2 the effect of changing the 
Haldane pseudopotentials was studied in Refs. [TBI [^. 
For a model with only the lowest two pseudopotentials, 
with ratio a — V2/V0, it was found that the groundstate 
is described by the bosonic Laughlin state for a < 0.4 
|23) . For a > 0.4 the groundstate is a compressible crys- 
talline phase related to the bubble crystal phase of mean- 
field theory described above. For larger filling factors, 
v > 1/2 the crystalline phase becomes increasingly sta- 
ble. Thus, for a > 0.4 the groundstates at all filling 
factors ly > 1/2 are compressible crystalline states. This 
has been confirmed in numerical studies for filling factors 
iy = l,3/2,2[2HH3n|. 

In the present model, with a non-Abelian gauge field, 
we have shown that for < /3 < v^ the lowest energy 
single particle state has n — 1 and only the pseudopo- 
tentials Vq and V2 are nonzero. The ratio a = V2/V0 is 
in the range < a < 0.118. Thus, over this range, one 
expects the groundstate at z^ = 1/2 to be an incompress- 
ible liquid that is well described by the Laughlin state j23j. 
Similarly, over this range, a is sufficiently small that the 
groundstates at ;/ = 1,3/2 and 2 are well described by 
the Moore-Read state, the k = 3 Read-Rezayi state, and 
the k = A Read-Rezayi state respectively [57] ■ Indeed, a 
small non-zero value of a has been shown to tune the 
system into a regime where these very interesting non- 
Abelian quantum Hall phases describe the groundstate 
accurately [IHHSni. 

For \/3 < /3 < a/S, the lowest energy state has n = 2. 
The non-zero pseudopotentials are Vo,V2 and V4. Over 
this range of /3, the ratio a = V2/V0 ~ 0.65, while 
Vi/V{) < 0.2. Neglecting the effect of this small value 
of V4, we can again make use of the results for the pure 
Vb-V2 model. (Neglecting V4 was accurate for interpret- 
ing the mean-field groundstate described above.) Now, 
the ratio a = 0.65 is so large (> 0.4) that the ground- 
state at all filling factors v > 1/2 is expected to be a 
compressible crystalline state |23lE51 - l5D] . 

Thus we can conclude that, for the filling factors 
v = 1/2,1,3/2,2, there is a transition from incompress- 
ible quantum liquid states (Laughlin and Read-Rezayi 
states) to a compressible crystalline state as the non- 
Abelian field is increased through /3 = a/S- This conclu- 
sion is in agreement with independent exact diagonaliza- 
tion calculations reported in Ref.[15| in those parameter 
regimes for which there is overlap. 



VI. CONCLUSIONS 

We have studied vortex phases in a model for charged 
particles in a non-Abelian gauge potential pertaining to 
a U{2) symmetry. Applying Gross-Pitaevskii mean field 
theory we have considered the effect of contact interac- 



tions between particles. These lead to the formation of 
triangular, square, and bubble crystal lattices for increas- 
ing values of the parameter for the non-Abelian term in 
the gauge potential. We calculated the Haldane pseu- 
dopotentials for the energy states of the non-Abelian 
model and find that they indicate effective interactions 
with non-zero range in the system. We find a general 
agreement with results on the effect of long-range in- 
teractions [TK [23] . We conclude that the effect of the 
non-Abelian gauge field on the properties of the lowest 
energy state is equivalent to the introduction of effective 
interactions with non-zero range. 

We suppose throughout that the interactions are small 
compared to Landau level spacing. However, at values for 
the parameter /3 ~ ^/In -\- 1, where the successive energy 
levels for non-Abelian model cross, our approximation is 
not valid, since Landau level mixing [31] is then expected 
even for small interactions. 
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These wavefunctions are orthonormal within the peri 
odic cell. Consider, for example, the wave functions ( Al I 
We have 



"'0 



'/'o,fci'/'o,fe2da;dy = 



E 

Pl,P2 



Ski+piNv,k2+P2Nv f " p-iYk+pLy-yf 
7rl/2 



(A4) 



where all summations in the p symbols (here and in the 
following) extend from — oo to oo. The S symbol in the 
latter equation gives zero for every pi ^ ^ p2 since the 
k's take values in the range < fci, /c2 < Ny. Therefore 
it is equal to the product Sk-^,k2 ^pi,p2- Using the result 



{Y^+pLy^yf _ / ^„^-iy~Y„,f 



E / '^y- 

„ "'0 



dye' 



'TT, 



which we substitute in Eq. ( A4 ) , we finally find that 00 



(A5) 



are orthonormal. A similar procedure for the first (n = 1) 
and second (n = 2) Landau level wave functions proves 
the orthonormality condition 



Appendix A: Spatially periodic wave functions 

The single-particle energy eigenstates of the non- 
Abelian gauge field are given by Eqn. ([7|, in which (j)n,k 
are the normalized Landau level wavefunctions for uni- 
form Abelian magnetic field. For the periodic system we 
study, with a rectangular cell of size L^ x Ly, these func- 
tions must have spatial periods of L^ and Ly in the x 
and y directions. 

The lowest Landau level wave functions with these pe- 
riodicities are ED 



PO.fc 



^ oo 

1/2U/2 2^ 



(i.7ri/2)i/ 



i{Yk+pLy)x -^{Yk+pLy~y)'^ 



-'yJ-^ g 2 



(Al) 



Yk ^ ^^K 



where k takes the integer values fc = 0, . . . , TVy — 1, with 
Ny given by Eqn. (11). Similarly, the spatially periodic 



wave functions in the first and second Landau levels are 



n,k 



(L. 2^1/2)1/ 



^ OO 

— -— E '^{Yk+pLy-y)e 



e 2 



\{Yk+pLy-yf 



i{Yk+pLy)x 



(A2) 



KM't^"'i'2 dydx = 5k^,k2 



(A6) 



Appendix B: Interaction potentials 

The interaction potentials entering in the calculation 



of the interaction energy ( 13 1 are of the form 



Vi, 



ki,k2,k3,ki 



where 



TT-l i -I i 

ki,k2,k3,k4 



V, 



t t t t — 



V, 



ki,k2,k3,k4 

T ; i T - 



^1,^2, ^3,^4 



V, 



I t t -I 



ki,k2,k3,ki 



V, 



■J' ■4' -4' ■4' 



ki,k2,k3,ki 



/?» + l 

2/?„ 

a 11^ ki,k2,k3,ki ' ki,k2,k3,ki' 

/3n-l\ T.t t t T 






(Bl) 



KM'^nM^"M(t>n,ki d^r, 



and we have explicitly taken into account that the outgo- 
ing particles (fci, fc2) have the same spin as the incoming 



ones (^3,^4). We present results for n = 1 and n = 2 in 
this paper. 
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